to appear in ApJ 

A First Glimpse into the far-IR properties of high-z UV-selected 
Galaxies: Herschel/PACS observations of z~3 LBGs 

O : G.E Magdis\ D.Elbaz^ H.S. Hwang ^ E. Daddi\ D. Rigopoulou^ B. Altieri^ P. 

Andreani^, H. Aussel^, S. Berta^, A. Cava^, A. Bongiovanni^, J. Cepa^, A. Cimatti^, M. 

Dickinson^, H. Dominguez^, N. Forster Schreiber^, R. Genzel^, J-S Huang ^^, D. Lutz^, R. 
Maiolino^^^ B. Magnelli^ G.E. Morrison ^°'^\ R. Nordon^ A.M. Perez Garcia^ A. 
^ ■ Poglitsch^, P. Popesso^, F. Pozzi^^, L. Riguccini^, G. Rodighiero^^, A. Saintonge^, P. 



(N 



o 
u 

6 



C^ 



Santini^^, M. Sanchez- Portal^, L. Shao^, E. Sturm ^, L. Tacconi^, I. Valtchanov 



ABSTRACT 



We present first insights into the far-IR properties for a sample of IRAC 
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c/5 ! and MIPS-24/im detected Lyman Break Galaxies (LBGs) at 2; ~ 3, as derived 

from observations in the northern field of the Great Observatories Origins Survey 
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(GOODS-N) carried out with the PACS instrument on board the Herschel Space 
Observatory. Although none of our galaxies are detected by Herschel, we employ 
a stacking technique to construct, for the first time, the average spectral energy 
distribution of infrared luminous LBGs from UV to radio wavelengths. We derive 
a median IR luminosity of Ljr = 1.6 x 10^^ Lq, placing the population in the class 
of ultra luminous infrared galaxies (ULIRGs). Complementing our study with 
existing multi-wavelength data, we put constraints on the dust temperature of 
the population and find that for their Ljr, MIPS-LBGs are warmer than submm- 
luminous galaxies while they fall in the locus of the LiR-Td relation of the local 
ULIRGs. This, along with estimates based on the average SED, explains the 
marginal detection of LBGs in current sub-mm surveys and suggests that these 
latter studies introduce a bias towards the detection of colder ULIRGs in the 
high-z universe, while missing high-z ULIRGS with warmer dust. 



1. Introduction 

Lyman break galaxies (LBGs), represent a significant fraction of z ~ 3 star- forming 
galaxies that have been discovered up to now. Since LBGs are the most common popula- 
tion at this epoch, investigating their properties is crucial to enhance our view of the early 
universe. A large number of mult i- wavelength studies, spanning from X-rays to radio, have 
characterized the nature of the population. One recent advance into our understanding of 
this population was provided by the Infrared Array Camera (IRAC), the Multi-band Imag- 
ing Photometer (MIPS) and the Infrared Spectrograph (IRS) on board the Spitzer Space 
Telescope (Spitzer). In particular, IRAC measurements at rest-frame near-infrared wave- 
lengths (3.6-8/im) indicate that their stellar masses are typically lO^-lO^^M© (e.g., Magdis 
et al. 2010a) while the detection of a fraction (~20%) of these UV-selected galaxies at MIPS- 
24/xm unveiled the sub-population of the infrared luminous Lyman break galaxies (Huang et 
al. 2005, Rigopoulou et al. 2006). Subsequent follow-up IRS spectroscopy of MIPS detected 
LBGs, revealed prominent Polycyclic Aromatic Hydrocarbon (PAH) features in their mid-IR 
spectrum, providing evidence for the existence of significant amounts of dust and that they 
are mainly powered by star-formation rather than an AGN (Huang et al. 2007, Rigopoulou 
et al. 2010). 

The star-formation rate (SFR) of MIPS detected LBGs has been studied in detail by 
Magdis et al. (2010b), where they found that UV, mid-IR and radio measurements of the SFR 
for UV-selected galaxies are in good agreement, indicating an average <SFR>=250MQ|/r^^. 
This SFR along with the fact that these galaxies contain considerable amounts of dust 
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as indicated by their mid-IR emission, suggests that MIPS-LBGs should be detectable at 
sub-mm wavelengths. Although the first attempts to detect far-IR counterparts were not 
successful (e.g. Chapman et al. 2000), recently there has been a significant progress in this 
front. Chapman et al. (2009) reported the SCUBA-850/im sub-mm detection of Westphal- 
MM8 while Rigopoulou et al. (2010) reported MAMBO 1.2mm detections of a further two 
LBGs, EGS-D49 and EGS-M28 selected based on their strong MIPS 24/im emission. In 
addition to those, two more lensed-LBGs have been detected in the sub-mm by Baker et al. 
(2001) (MS1512-cB58) and by Siana et al. (2009) (Cosmic Eye). Despite these promising 
detections though, the FIR properties of LBGs still remain elusive. 

In this paper we use observations of the northern field of the Great Observatories Origins 
Survey (GOODS-N) obtained by Photoconductor Array Camera and Spectrometer (PACS) 
(Poglitsch et al. 2010) on board the Herschel Space Observatory (Pilbratt et al. 2010). We 
complement our study with multi-wavelength data and we aim to get a first insight into the 
FIR properties of a sample of MIPS 24/im detected LBGs, construct the average spectral 
energy distribution (SED) of the population, derive estimates of their dust temperature (Td) 
and compare their properties to those of SMGs. The Herschel data used in this study are 
part of the Herschel Science Demonstration Phase of the PACS Evolutionary Probe (PEP, 
PI Lutz) program. Throughout this paper we assume r2jn=0.3, Ho=71km sec~^ Mpc~^ and 
fiA=0.7 



2. Sample Selection and Herschel Observations 

GOODS-North was observed by PACS, at 100- and 160/im in the framework of the 
Guaranteed Time Key Program, Pacs Evolutionary Probe (PEP). The area covered by PACS 
was lO'xlS' in a total exposure time of 30 hours and an average integration time of ~ 2.5sec 
per pixel. For source extraction two independent methods were employed: the traditional 
blind source extraction that was performed by using the Starfinder Point-spread function 
(PSF)-fitting code (Diolaiti et al. 2000) and a guided extraction using 24/im priors. In 
both methods, PSFs were extracted from the final science maps and were subsequently used 
for a PSF-fitting photometric analysis. Flux uncertainties were obtained by Monte-Carlo 
simulations, adding fake sources in the real maps. The measured la noise is 1.00 mJy at 
100/im and 1.90 mJy at 160/im. For more details we refer the reader to the appendix of 
Berta et al. (2010). 

Our sample of LBGs is identical to the one studied by Magdis et al. (2010b) and we 
refer the reader to this paper for a detailed description. In brief we consider 49 IRAC de- 
tected ([3.6]ab <25.0 or 0.36/iJy) LBGs in GOODS-N, that have originally been optically 
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selected (U,G,R, R<25.5) by Steidel et al. (2003). All galaxies in our sample have spectro- 
scopic redshift (median z=2.95) derived by optical spectroscopy that has also been used to 
determine the absence of AGN signatures (i.e., strong high ionization emission lines) in their 
rest-frame UV spectrum (Steidel 2003, Shapley 2003). Among the original sample, 9 LBGs 
are also detected at MIPS-24/im {5a ~20/xJy) with a median 24/im flux density /24=31.2/iJy. 
Henceforth, we will refer to these LBGs as the MIPS-LBGs while to the remaining (40) as 
IRAC-LBGs. IRAC photometry for these objects is provided in Magdis et al. (2010a), 
while an analysis of MIPS 24/im data is described Magdis et al. (2010b). Matching our 
sample with the 100- and IGOfim. blind and prior catalogs returned no individual detections 
at the 3 sigma level (3mJy and ~5mJy respectively). Since none of our LBGs is individu- 
ally detected we use the residual maps to do stacking. On stacking we first considered the 
sample of MIPS-LBGs. We employed median stacking analysis, cutting sub-images of the 
residual maps centered at the optical position of each undetected LBG. The residual maps 
were created by Starfinder, by subtracting all individually detected sources down to 3a. To 
avoid contaminating the stacked signal from residuals, we only added galaxies to the stack if 
there were no bright sources within ~4"(100/im) and ~6"(160/im). Then a stacked flux was 
measured in a manner similar to the measurement of the detected PACS sources. To quan- 
tify the error of our measurement we stacked at 9 random positions and repeated it 50.000 
times. The la of the distribution of the derived fluxes was adopted as the uncertainty of our 
measurement. Stacking at 100/im returned no detection (fioo=0.47±0.36 mJy) indicating a 
3a upper of fioo <l.lmJy, while at 160/im we recovered a median flux density /i6o=2.21 ± 
0.52mJy (S/N~4.2) for the MIPS-LBGs. The final 160/im stacked image is shown in Figure 
1 along with the distribution of the fluxes that correspond to random stacked positions. 
This distribution indicates the probability to recover a stacked flux density /leo > 4.2mJy 
by chance, is 2.31 x 10^^. Stacking the whole sample of LBGs, or just the IRAC LBGs alone, 
returned no detection in either PACS band. For the IRAC LBGs we determine a 3a upper 
limit of fioo < 0.47mJy and figo < 0.79mJy while the actual fluxes that correspond to the 
stacked images are /ioo=0.21±0.16 mJy and fi6o=0.32± 0.26 mJy. 

For the IRAC sample we also consider median stacked flux densities of fi.imm=0.41±0.11mJy 
and fi.4GHz=3.6±0.8//Jy as derived from stacking at 1.1mm AzTEC and 1.4GIIz VLA map, 
while for the MIPS sample we consider a median stacked flux density of fi.4GHz=8-5±2.2/iJy 
and a 3a upper hmit of fi.imm=l-lniJy (Magdis et al. 2010b). We also add optical (ground 
based UGR, ACS BViz) and near-IR (J,K) photometry provided by the GOODS team. 
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3. Results and Discussion 

3.1. Spectral Energy Distribution of LBGs 

Using the mult i- wavelength photometry (and the upper limits) described above we at- 
tempt to construct the full, (rest-frame UV to radio) average SED of MIPS selected LBGs. 
We fit the optical to near-lR part with model SEDs generated by the Chariot and Bruzual 
2007 (CB07) code and the mid-IR to radio with template SEDs from the Chary & Elbaz 
2001 (CEOl) and Dale & Helou 2002 (DH02) libraries. For the fitting procedure we allow 
for renormalization of the CEOl templates while for the DH02 models the luminosities are 
normalized as described by Marcillac et al. (2006). Results based on the two methods are 
in very close agreement indicating a median Lir=1.6(±O.5)x1O^^L0 and that MIPS-LBGs 
belong to the class of Ultra-Luminous Infrared Galaxies (ULIRGs). This luminosity trans- 
lates to SFR= 275MQyr-^ (Salpeter IMF), in good agreement with the UV (25O+^^M0yr-i) 
and radio (28O±85M0yr^^) SFR estimates presented in Magdis et al. (2010b). Adopting 
a median stellar mass of MIPS LBGs M*=7.9xlO^°M0 (Magdis et al. 2010a), we derive a 
specific star-formation rate of SSFR=3.5Gyr~^. This value is lower than that derived for 
the IRAC-LBGs (4.3Gyr~\ Magdis et al. 2010a) but higher than that found for z~2 star- 
forming galaxies (2.5 Gyr~^, Daddi et al. 2007a). Finally, based on the best fit CEOl model, 
we estimate a median /85o=1.36mJy, very close to the confusion limit of current sub- mm 
surveys. We comment on this finding later on where we present a comparison between LBGs 
and SMGs. 



3.2. Indication of Warmer Dust in MIPS-LBGs 

As we have already discussed, stacking IRAC-LBGs at 1.1mm has recovered a fiux 
density of /i.imm=0.41±0.1mJy while at 160/im there is no detection. The situation is exactly 
the opposite for the MIPS-LBGs where stacking at 1.1mm indicate no detection while at 
160/im we recover a median fiux of /i6o=2.21±0.52 mJy. This inversion could be interpreted 
as a change in the shape of the SED of the two samples. Furthermore a detailed study of 
IRAC-LBGs by Magdis et al. (2010b), indicated that the average IRAC detected LBG is 
a Luminous Infrared Galaxy (LIRG) with Lir=4.5(±1.2)x10^^ Lq, less luminous than the 
current sample of MIPS-LBGs. It has been observed that in the local universe ULIRGs have 
warmer dust and their SEDs peak at shorter wavelengths than those of LIRGs. Here we 
attempt to put constraints on the peak of the SED and the dust temperature (T^) of LBGs 
and investigate these parameters as a function of luminosity at z~3. 

Although ideally such a task requires photometric points close to the peak of the SED 
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of a galaxy (which in this case translates into SPIRE data), here we use PACS 160/im and 
AzTEC 1.1mm flux densities and upper limits, with the advantage of sampling both sides 
around the peak of the SED. Using 3a upper limits for the /igo of IRAC-LBGs could provide 
us with an upper limit of the T^ of the sample while 3a upper limits for the /i.imm of the 
MIPS LBGs will yield an estimate of a lower limit of the T^ of the sample. 

We first employed Monte Carlo simulations to generate 10.000 combinations of 160/xm 
and 1.1mm fluxes. For each of the two bands, the simulated fluxes were generated following 
the flux distribution derived from the stacking simulations and centered on the measured 
stacked flux. To derive the dust temperature for each of the generated realizations, we used 
a single temperature greybody fltting form, Fiy oc u^'^^ / {expihu / kTd) — 1), where /3 is the 
dust emissivity and T^ is the effective or emission weighted dust temperature. We flxed the 
P value to 1.5 (Gordon et al. 2010) and considered the median redshifts as representative 
for each sample (z=2.98 for IRAC-LBGs and z=2.92 for the MIPS-LBGs). The T^ for 
each realization was then obtained from the best flt model that was derived based on the 
minimization of the x^ value. Since at lOO/im (rest-frame 25/im) there is strong contribution 
of Very Small Grains (VSGs) in the IR emission, /loo upper limits were not fltted but they 
were considered as a sanity check, i.e. that the models do not violate the upper limits. 
Based on the distribution of the obtained dust temperatures (Figure 3) we then derive an 
upper limit of T^ = 38.2(39.3) K (IRAC-LBGs), and a lower limit of T^ = 40.1(38.9) K 
(MIPS-LBGs) at a 2a{3a) confldence level. The best flt models that correspond to the 2a 
temperature limits are depicted in Figure 4. We note that adopting /3 = 2.0 instead of 1.5, 
would result in a lower dust temperature, by 2-3 K. 

The above analysis indicates that the dust temperatures of the two samples are different 
at a 2.6a confldence level, with MIPS-LBGs having warmer dust. Since MIPS-LBGs and 
IRAC-LBGs are objects divided by a threshold in their rest 6/im emission, which is due to 
warm, very small grains and/or PAHs, if the presence of warm small grains correlates with 
the presence of warm larger grains then our result would be object of selection bias. To ex- 
plore this possibility, we consider a sample of local ULIRGs with ISOCAM 6.75/im and IRAS 
60- and lOO/xm data (Elbaz et al. 2002). Plotting the /loo/Zeo which is a dust temperature 
indicator, against the rest-frame 6.75/im luminosity (that corresponds to 24/im emission of 
z~3), reveals that there is no obvious correlation between 6//m luminosity and Tj. Hence, 
assuming that this applies for galaxies at z~3, it seems that our result is not subject to a 
selection bias. On the other hand, as we have already shown, MIPS-LBGs are more infrared- 
luminous than IRAC-LBGs with corresponding infrared luminosities Lir=1.6(±0.5) xIO^^Lq 
for MIPS- and Lir=4.5(±0.5)x10"Lq for IRAC-LBGs. Therefore, it appears that we ob- 
serve for the flrst time the general Lir— Td trend seen in the local universe, (with more 
luminous galaxies having warmer dust), for UV-selected galaxies at z ~ 3. 
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3.3. Comparison with SMGs 

Considering tlie large SFRs and substantial dust reddening that are inferred for some 
LBGs, and the fact that their spectra exhibit characteristics of local starbursts (Pettini et 
al. 2001), it is somewhat surprising that there are only few examples of direct sub-millimeter 
detection for these galaxies. MIPS-LBGs are the most rapidly star-forming, most luminous, 
and dustiest galaxies among the high redshift UV-selected population, and therefore are 
the best candidates for having far-IR emission that could be detected in current sub-mm 
surveys. In what follows we discuss how the current data can shed light on the mystery of 
this marginal detection of LBGs in sub-mm bands. Based on the average SED of MIPS- 
LBGs that we presented here, we predict that the flux density of the MIPS-LBGs emitted 
at 850/im is /85o=l-l-l-5mJy, just below the current confusion limit. It could therefore be 
suggested that MIPS-LBGs provide a link between SMGs and typical UV selected LBGs 
that are faint in the IR. It has also been speculated that the fact that LBGs are not detected 
in the sub-mm bands is due to warmer dust compared to that of SMGs (Kaviani et al. 2003). 
Although to test this one needs both robust measurements of the Tj (based on observations 
that probe the peak of the SED of the galaxies) and to have a large sample of LBGs, the 
upper limits that we derived in this study can provide a first insight. 

Figure 5 compares dust temperature versus infrared luminosity for the MIPS and IRAC 
LBGs as derived from this study. For comparison we also consider the large compilation 
of z ~ 2 SMGs by Chapman et al. 2005, and the local/intermediate-z samples of ULIRGs 
presented by Clements et al. (2010) and Yang et al. (2007). In all these studies, the T^ was 
measured in manner that is similar to ours, fitting modified black-body models with fixed 
/9=1.5 to the FIR photometric points of each sample, hence the comparison between these 
studies and our sample is meaningful. We also plot, the 3a envelope of the Li^i-T^ relation 
for local infrared galaxies in SDSS (Hwang et al. 2010 in prep). 

It is evident that for the Ljr of the MIPS-LBGs, the bulk of SMGs are considerably 
colder, while MIPS-LBGs fall in the locus of the local ULIRGs and are within the scatter 
observed in local galaxies. Based on modified black body models, we also compute tracks 
of constant 850/im flux density for galaxies at z=3, close to the confusion/detection limit of 
current sub-mm surveys (/85o=lniJy and /85o=2mJy). MIPS-LBGs lie in between the two 
tracks, indicating that a typical MIPS detected LBG emits at 1-2 mJy level at the sub-mm 
bands, in excellent agreement with the median /gso flux density derived from the average 
SED of the population. This explains the small overlap between the LBGs and SMGs found 
in previous studies. For a given star formation rate and optical depth for dust absorption 
(and hence, for a given bolometric far- infrared luminosity), the warmer dust temperature 
in LBGs leads to fainter sub-mm emission compared to that found in typical SMGs. This 



finding also suggests tliat star-formation in LBGs is spatially compact. On the other hand, 
recent studies based on minimal excitation (CO J=l -^ 0) emission that traces the bulk 
of the metal-rich molecular gas provide evidence of a more extended distribution of star- 
formation in SMGs (e.g. Ivison et al. 2010), that could explain their colder dust component. 
Generalizing our results, we argue that ULIRGs with warmer dust than that of SMGs exist 
in the high-z universe and would not be selected in current sub-mm surveys. This is further 
enforced by the study of IRAC selected ULIRGs at 2; ~ 2, (Magdis et al. 2010 in prep), 
where based on Herschel PACS/SPIRE observations, we find that 60% of the sample would 
be missed by current sub-mm surveys. 

We have presented a first glimpse into the far-IR properties of z ~ 3 LBGs, demonstrat- 
ing the power of Herschel in the study even of high-z UV-selected galaxies and its ability to 
address issues that we were unable to resolve in the pre-Herschel era. We aim to further ex- 
tent this study, complementing PACS with SPIRE data that probe the peak of the SED of z 
~ 3 galaxies and taking advantage of the Super-Deep and Ultra-Deep surveys in GOODS-N 
and GOODS-S respectively as part of the Herschel-GOODS program (PI D. Elbaz). 

PACS has been developed by a consortium of institutes led by MPE (Germany) and 
including UVIE (Austria); KUL, CSL, IMEC (Belgium); CEA, OAMP (France); MPIA 
(Germany); IFSI, OAP/AOT, OAA/CAISMI, LENS, SISSA (Italy); lAC (Spain). This 
development has been supported by the funding agencies BMVIT (Austria), ESA-PRODEX 
(Belgium), CEA/CNES (France), DLR (Germany), ASI (Italy), and CICT/MCT (Spain). 
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Fig. 1.— Left) PACS-160/im stacked image (130" x 130") of 9 MIPS-24/im detected LBGs in 
GOODS-N. Right) Stacking simulations at 160/im. Distribution of the measured fluxes 
derived from 50000 stackings at 9 random positions along with the best Gaussian fit 
(rms=0.52mJy). The red line denotes the flux measurement of the stacking at the posi- 
tion MIPS-LBGs indicating a ~4.2cr detection. 
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Fig. 2. — Rest-frame average SED of z~3, MIPS-LBGs. For the SED we use the median 
UGR+BViJK+IRAC+MIPS24 photometry of MIPS-LBGs, and the values (or upper hmits) 
derived from stacking 100/xm, 160/im AzTEC and radio. Upper hmits are indicated with 
green squares. The rest-frame UV-NIR portion of the data is overlaid with the best fit CB07 
model (orange line), while the mid-IR to radio is shown with the best-fit CEOl model (blue 
line). Black diamonds show the fluxes predicted from the best fit CEOl model. 
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Fig. 3. — Distribution of the derived T^ for the IRAC- and MIPS-LBGs based on MC simula- 
tions. The two distributions begin to overlap at a 2.6a confidence level. Blue lines illustrate 
a gaussian fit to the left (MIPS-LBGs) and right (IRAC-LBGs) part of the distribution. 
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Fig. 4. — Derivation of lower and upper T^ limit estimates for MIPS- and IRAC-LBGs 
respectively. Filled black squares indicate the median stacked flux density of MIPS-LBGs 
at 160/im and the 3a upper limits of the flux density at 1.1mm as estimated from stacking. 
The solid red line shows the best fit modified black body model that corresponds to the 2a 
Td lower limit. For the model we assume fixed /3=1.5. Similarly, filled green circle indicate 
the median flux density of IRAC-LBGs at 1.1mm from Magdis et al. 2010b and the 3a flux 
upper limit at 160/im as derived from stacking in the present study. The solid blue line 
shows the best flt modifled black body model that corresponds to the 2a T^ upper limit. 
The open black squares and green circle indicate 3a flux upper limit at 100/im as derived 
from stacking for MIPS- and IRAC-LBGs respectively. Vertical black dotted and dashed- 
dotted lines indicate the wavelength where the far-IR emission peaks for MIPS-LBGs and 
IRAC-LBGS. 
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Fig. 5. — The Lir— T^ relation for MIPS and IRAC LBGs (red and black circles respectively). 
We note that these values correspond to 3a upper and lower T^ limits for IRAC and MIPS 
LBGs. Included are results for 0<z<0.98 ULIRGs (blue filled triangles, Clements et al. 2010, 
Yang et al. 2007). The green shaded region depicts the loci of high-z SMGs by Chapman 
et al. (2005), while the orange shaded area shows the 2a envelope of the Lm — Tj relation 
for local IR galaxies in SDSS (excluding AGNs) adopted from Hwang et al. (2010 in prep). 
Black dashed lines represent tracks of constant flux density at 850/im (/85o=lmJy and /85o=2 
mJy) for galaxies at z=3. Objects at z=3 with higher /gso, lie on the right of the lines. This 
plot demonstrates that for a given Ljr, the bulk of SMGs are considerably colder than MIPS- 
LBGs which fall in the locus of the local ULIRGs. We also note that for their T^ and Lir, 
MIPS-LBGs fall in between the 1- and 2 mJy constant flux density at 850/im, indicating 
that they would be missed by current ground-based sub-mm surveys. Finally, there is a hint 
that the LBGs at z ~3 follow the general trend observed in the local universe, with more 
luminous galaxies having warmer dust. 



